ADA104877 


E0ARD*TR-81  -  1  Q 

'k 


SCHOOL  Of  ENGINEERING 
TEL-AVIV  UNIVERSITY 


Department  of  Eletronlc  Communications, 
Control  and  Computer  Systems 
Ramet-AvW,  Tel-Avlv  69978,  Israel 


DISTRIBUTION  STATEMENT  A 

v  - 

Approved  lor  public  release} 

Distribution  Unlimited _ 

819  2  8 

% 


DTIC 

ELECTE 
SEP  30  1981 

D 


2  1 8 


EOARD-TR -8  1  -  1  0 

Report  Numb i  \^FOSK-80-0 178^002 


_Acri5Si0n  For 

NT  IS  GRA&I  sV" 

DTIC  TAB  V,  I 

Unannounced  n  | 

Justification _  I 


By _ _ _ 

Distribution/ 
Aviiilabilitv  Curb's 
.Avail  iiiid/or 
)ist  j  Special 


FTXKD  GAIN  CONTROLLFR  ni'.STGN  FOR  ATRGRAFT. 


Ynrnm  Baram  and  Duron  '  V.  i  de  Iman 


department  of  Fleet roni .  Systems 
Sell i .o  1  oi  F.ngincer Iny. 


Tel-Aviv  University 


//  AOtb  dun  *-,49  81 


I,- 


Final  Scientific  Re'p«rt,’  10  Jun<d98$  —  31.0  May  W81,  (  ci  en  t  I  f  i  < 


Approval  lor  public  release,  distribution  unlimited 


-  v! 


Prepared  for  : 


U.S.  Air  Force  Office  of  Scientific  Rescan 


Furopean  Office  of  Aerospace  Research  and  Development,  London,  F.ngland. 


'i'/Q. 


RE  iK'Ci*11'1*. N  !  r\  1  lO.t 


j  IIK t’ORT  l»OC» V'KN' VAT iO!i 

READ  1 :  ’.'>  i':;  HGTION3 

f-ewpe  c.::.”  rrw;  rosi: 

:  J  .  .  i  Mu: Tier 

1  cU..:',i!-]:i-R  |-10 
i 

2.  (P'vv  rre.-.u:.  !  .--i  No. 

i.  S'  .  *  jiielif.  Cft  a  log  ’e.-ilier 

j  t<  t  );■•  j.ir:  Sue  i' tie  ) 

M-.s-i  Gain  Control  irr  Css i on  Cor  Aircraft 

5.  Type  of  Report  A  Period  Covered 

final  Scientific  Report 

1.  Juj.e  1'jhO  -  31  ,.'ijy  19!!  1 

b  .  Per  forming  C:  .  Report  Number 

r  STR  SI  -Of) 

J  7  .  Author ( s ) 

!  Dr.  Yoram  Baram  and  Doron  EMclnwn 

i 

.1.  'out rant  or  Grant  Number 

Gr.mt:  A  RISK  -RO-0178  .  (L 

'J  Dor  forming  Organisation  Sane  and  Address 

Department  of  electronic  Systems 

School  of  hop  inner ing 

Te 1 -Av i v -Un i vers i ty 

10.  program  Element,  Project ,  Tank 
Area  6  Cork  Ur  It  Numbers 

6 .110  2  F 

2301-D1 

\ 

11.  Controlling  Office  Name  and  Address 

European  Office  of  Aerospace  Research 
!  and  Development 

i  Box  14,  FPO  New  York  09510 

12.  Report  Date  (1 

30th  Juno  ]%1.  |j 

......  3 

12.  Number  of  Pages  j 

^  1 

j  !  C  Monitoring  Agcncv  Name  and  .V  dress 

1 

1 

! 

ii 

I 

id.  A  1/  Distribution  Statement 

i 

Approved  lor  public  release;  distribution  unlimited.  r 

i  id  Supp  1  iwncary  Notes  I 


aircraft  control,  adaptive  control,  fixed  gain  controllers,  control 
systekn  filters 


.'ll  I: 


A  Method  for  de«  ic;rin«*  fixcJ  qe  •:  o  controllers  and  filter.,  for 
system-,  with  largo  paro:?oLo>-  vjt lotion  is  presented. 

The*  approach,  based  of .r.nr.PVM  ir foliation  criteria,  is  used  to 
design  a  non-adapti ve  bach-up  control  systen  for  a  given  aircraft 
and  is  shown  to  provide  good  performance  qualities. 


ccARo-TR-ei 

This  report  has  been  reviewed  by  the  FOARD  Information  Office  and  is 
releasable  to  the  National  Technical  information  Service  (NTIS).  At. 

N"'J.S  it  will  be  releasable  to  the  general  public,  including  foreign 
nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


✓ 

.// 


/  /'  //  )  /? 

t  ".JOHN  W.  DAfLET 
Major,  USAF 
Chief,  Electronics 


FOR  THE  COMMANDER 


CORDON  L.  HERMANN 
it  Colonel,  USAF 
Deputy  Commander 


A.%.. 


/ 


Introdue  t ion 


A  major  difficulty  associated  with  aircraft  control  is  that  the  dynamic 
characteristics  are  highly  dependent  on  variable  parameters,  such  as  the 
altitude  and  the  Mach  number.  Recent  years  have  seen  considerable  effort 
towards  the  design  of  adaptive  control  systems,  taking  the  dynamic  dependence 
on  varying  flight  parameters  into  account.  Yet,  the  design  of  simple  controllers 
and  filters,  minimizing  the  necessary  airborn  computation  and  data  storage 
capabilities,  on  the  one  hand,  and  the  dependence  on  accurate  air  data,  on  the 
other,  remains  of  major  interest.  Fixed  gain  controllers  and  filters  are  commonly 
used  as  back-up  control  systems  for  advanced  high  performance  aircraft. 

In  this  paper  a  method  for  designing  fixed-gain  controllers  and  estimators 
for  systems  with  large  parameter  variations  or  uncertainties  is  described.  The 
approach  is  based  on  minmax  criteria  defined  on  Rollback's  information  measure  [ 1 1 
and  previously  used  for  solving  model  simplification  problems  ([2],  [3]).  (We 
note  here  that  other  minmax  design  criteria  for  deterministic  systems  have  been 
proposed,  e.g.  [a]).  The  design  objective  is  that  the  maximum  possible  difference 
between  the  optimal  system  at  the  actual  operating  point  (i.e.  the  optimal  adaptive 
system)  and  the  selected  fixed-gain  system  be  minimal.  Finite  time  and  regulator 
design  criteria  are  defined.  The  proposed  technique  is  used  to  design  a  longitudinal 
back-up  control  system  for  a  given  high  performance  aircraft.  Simulated  aircraft 
responses  and  a  pilot  rating  chart  indicate  good  performance  qualities  of  the  resulting 
fixed  gain  control  system  for  the  entire  flight  regime. 
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2 •  Design  Criteria 


A.  General 


In  the  design  problems  under  consideration  we  have  a  stochastic  system 
M(p,s)  depending  on  a  parameter  p,  whose  values  may  vary  over  a  set  P  and  a 
parameter  s,  whose  value  is  to  be  selected  from  a  set  S  in  the  design.  Employing 
some  optimal  design  criterion,  to  each  parameter  value  p  there  corresponds  an 
optimal  value  s(p)  of  s.  It  is  desired,  however,  to  find  a  single  value  s° 
that  will  be  optimal  in  some  sense  for  the  entire  parameter  set  P.  The  cost 
(or  damage)  of  the  resulting  design  at  a  given  parameter  value  p  is  defined  as 
the  distance,  in  some  measure,  between  M(p,s)  and  M(p,s(p)).  Employing  Kullback's 
information  as  a  distance  measure  ([ 1]— [3J )  the  maximal  cost  for  a  given  choice  of  s 
is 


max 

peP 


ES<P,loS 


fS(p) (Y) 
-JL _ 

fS(Y) 

P 


where  Y  is  a  set  of  observations  fS^(Y)  and  fS(Y)  are  the  probability  densities 

'  p  p 

s  (d  ) 

of  Y  corresponding  to  M(p,s(p))  and  M(p,s),  respectively  and  p  denotes 
expectation  with  respect  to  M(p,s(p)).  The  design  objective  is  to  minimize  the 
maximal  cost,  i.e. 


s(  }  fS(p)(Y) 

min  max  E‘  p  log  — ^ - 

seS  pcP  P  fS  (Y) 


/ 

/ 
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B.  Controller  Design  for  a  Linear  System  in  Finite  Time 
Consider  a  dynamical  system  in  the  form 


x 


n+1 


=  F  (p,s)x  +  B  (p,s)ti  +  G  (p)w 
n  n  n  n  n  n 


(2.1) 


'n  =»  H  (p)x  +  v 
n  n  n 


whore  x  is  Gaussian  with  moan  m  (o) 
o  o 

(non-feedback)  deterministic  input, 
sequences  with 


and  covariance  1’  (p)  ,  u  is  an 

o’  a 


w  and  v  are  zero-mean  Gaussian  white 
n  n 


external 

noise 


cov  (w  )  =  Q  (p) 
n  n 

cov  (v  )  =  R  (p) 
n  n 

p  is  a  variable  parameter  vector  taking  values  in  a  set  P  and  s  is  a  vector  of 

controller  gains  taking  values  in  a  set  S.  Any  linear  state  feedback  is  accounted 

for  in  the  matrix  F  (p,s).  We  denote  the  above  svstem  by  M(p,s). 
n 

Employing  some  optimal  control  criterion, to  each  parameter  value  p  there 
corresponds  an  optimal  gain  value  s(p)  and  the  corresponding  system  M(p,s(p)). 

Tt  is  desired,  however,  to  find  a  constant  gain  s°  that  will  be  optimal  in  some  sense 
for  the  entire  parameter  set  P.  As  in  the  case  of  optimal  controller  design  at  a 
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slven  parameter  value,  the  design  may  be  aimed  at  optimising  the  system's  response 
to  some  input  seguence  or  tc  an  initial  condition  In  finite  time.  let  us  denote 


UP..)  -  I  <lP> 


r'o’cv", 


“  t  JUS  e  N 

N  p  fs  (Y  ) 


N  N 

where  Y‘  =  (>’,,•••  >YN' 


For  the  Gaussian  case  under  consideration  we  have 


N  ! l  (p,s)l  i 

I  (P,s)  =  i  I  [log  — - - ~-~+  tr  Zn  (P*s)rn(p’s)'k] 

N  N  n=i  U  (p.s(p)); 


where  k  is  the  observations'  dimension  and 


'V.‘>  •  *1  ‘l-V>'n!n-l<>''S>,lVJ,n|n-l(p’S)1  > 


(p,„p»  -  ^(p>'ivn-y„1„.1iP'»<p»l‘V^|n-i'r-3(p)>1  ' 


.*>  ’  Ep<P>''lVynin-l<P",)UV‘n|n-l(!’'‘)1  ' 


V  '  ,  (p 

n ;  n- 1 


„s;  =  E;;yn!vn-b 


V  ,  ,(p.s(p))  =  E^(p){y  |Yn-1} 

n  n-1  P  n 
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For  the  given  system  Z  (p,s)  is  obtained  from 

n 


2:  (p,s)  =  H  (p)i>  (p,s)H  (p)  +  R  (p) 
n  n  n  n  n 


where  i  (p,s)  is  obtained  from  the  Ricatti  equation 
n 


2  (p,s)  =  F  (p,s)y  (p,s)FT(p,s)  +  G  (p)Q  (p)GT(p) 

n+1  rt  n  n  n  n  n 


-F  (p,s)K  (p  ,  s )  F  (p , s)K  T(p,s)FT(p,s) 
n  n  n  n  n 


R  (p,s)  =  C  (p,s)HT(p)F  (p » s ) 
n  n  n  n 


(Z  (p,s(p))  is  obtained  in  a  similar  manner,  replacing  s  by  s(p)).  V  (p,s)  is 
n  n 


obtained  from 


r  (p,s)  =  H  (p)U  (p,s)  +  X  (p  ,s)X  (p ,  s)  ]  H  (p)  +  R  (p) 
n  nn  n  n  n  n 


where  (p,s)  is  obtained  from 
ti 


<Mp.s)  =  F  (p,s)t  (p  >s)F  (p  ,s) 1  +  G  (p,s)Q  (p)  G  T(p,s) 
n  n  n  n  n  n  n 


in i t ia 1  ized  at 


(p)  0 

o 


i>  (p,s)  = 

o 


and  X  (p,s)  is  obtained  from 
n 


y  ( p , s ) X  (p,s)  +  B  (p,s)u 
n  n  n  n 


initial  i  .ted  it 


m  (p) 
o 


I  m  (p) 
L  o 


We  have  used  the  notat ion 


F  ( p  ,  s 
n 


fn(P)S) 


s  ( p  )  ) 


j  F  (p , s ) K  (p,s)H(p) 
n  rt 


F  ( p ,  s ) [ I -K  (p,s)H  (p)] 
n  n  n 


i;  •  .  ?_*  Jji"- 


Suppose*  chat  the  systca  is  Mrir  invariant,  L.e. 


:<  ,  •  l’(n,s)x  +  Rtp,s)u  +  (’(.n)w 

n-r  i  '  n  n 


y  ’  H(p)x  +  v 


-■  v  (v  i  *  |(p),  fn’MV  )  -  il<p) 


:.’»»*  external  •' ri.'s!- f  <>r  cl!'  h  ,  is  now  ;lf's'.ir''f|  const  a  a  I 


V TV  oil  r.1,1 

no  a  1  torn 

i;»l  iv* 

ill.  1  L\1 

St 

fr 

t  l»r  !  r’ 

'Mint  .it  i 

'  •  r  \ 

finite  t Into  design 

in  tills  e.i 

i  so 

stc-ifiv 

st; 

ite  oi 

f  retail 

vtOt  l!« 

•  s  \  y,  n 

result  Is  fundament 

.)!  to  I";.' 

*  j  Xt’C 

\  yti\  i.n  i  t- 

pte 

1.  n  1 1  u 

d  o s  i  J*  1 » 

mrt  iiod  . 

i 

!  Theorem  1 


Suppose  that  Cor  each  pep  and  each  si  S  the  system  (Fip,s),C(p),  l'(p))  is 
I'ont  rol  table  and  observable  and  the  controller  s(p)  Is  stabilizing  (i.e.  F  (p ,  p ) ) 
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has  all  its1  eigenvalues  Inside  the  unit  circle).  Then  the  cost  function 
l^tp.s)  has  a  finite  limit  value,  given  by 


I(p,s)  -  log  +  tr  r*l(p,e)r(p,e)-k  (2.6) 

I  £ (p.e(p) | 

where  I(p,s),  Z(p,s(p))  and  r(p,s)  are  the  limit  values  of  E  (p,s),  I  (p,s(p))  and 

n  n 

f  (p,s),  respectively, 
n 

Proof 

It  is  well  known  that  when  the  system  is  controllable  by  the  input  noise 

and  observable,  the  state  estimation  error  covariance  has  a  non-singular  limit 

value.  It  follows  that  the  output  estimation  error  covariance  £  (p,a)  has  a 

n 

non-singular  limit  value  E(p,s)  if  H(p)  has  independent  rows,  which  is  the  normal 

case.  It  is  also  well  known  that  an  equation  of  the  type  (2.4)  (a  Lyapunov  equation) 

has  a  limit  solution  for  1*  (p,s)  if  F(p,s)  has  all  its'  eigenvalues  inside  the  unit 

n 

circle.  But  due  to  the  structure  of  F(p,s),  its*  eigenvslues  are  the  eigenvalues 
of  F(p,s(p))  together  with  the  eigenvslues  of  F(p,s) [ I-K(p,s)H(p) ] .  Since  s(p>  i« 
stabilizing,  F(p,s(p))  has  all  its'  eigenvalues  inside  the  unit  circle.  Since  the 
system  (F(p,s) ,G(p) ,H(p)  is  controllable  and  observable,  the  matrix  F(p,s) [I-K(p,B)H(p) ] 
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ha  s  nil  its  eigenvalues  inside  the  unit  circle  (see,  e.g.  [5]  p.  144).  It  follows 

that  all  the  eigenvalues  of  h'(p,s)  are  inside  the  unit  circle  and,  consequently, 

(2.4)  has  a  limit  solution  for  i)1  .  In  addition  since  the  system  (2.5)  is  stable, 

n 

X  (p,s)  has  a  finite  limit  value.  It  follows  that  f  (p,s)  in  (2.3)  has  a  finite 
n  n 

limit  value,  r(p,s).  We  have  thus  shown  that  Z  (p,s)  and  T  (p,s)  have  limit 

n  n 

valuec  2(p,s)  and  r(p,s).  The  assertion  readily  follows  from  (2.2). 

It  should  be  noted  that  the  controllability  and  observability  conditions  are 
structural  requirements  which  are  "generically"  independent  of  the  specific  values 
of  the  parameter  p  and  s.  The  condition  that  s(p)  is  stabilizing  is  normally 
satisfied  in  pointwise  optimal  design.  For  instance,  consider  the  system 


■;  =  F  (p ) x  +D(p)z  +B(p)u  + 

n+1  n  n 


G(p)w 


=  H(p)x  +  vn 
n  n  11 


where  z  is  a  linear  state  feedback  process,  i.e. 
n 


(2.7) 


z  =  -sx 
n  n 


(2.8) 


Then  it  is  well  known  that  by  an  appropriate  choice  of  the  gain  matrix  s,  the 
eigenvalues  of  the  matrix 


F(s,p)  =  F (p )  -  D(p) p 


(2.9) 
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can  be  placed  inside  the  unit  circle  provided  that  the  pair  (F(p),  D(p))  is 
stabil  izable .  One  such  choice  of  s  is  obtained  by  the  linear  quadratic 
cri terion 

go 

min  [x  ^Lx  +  x  7.  +  z  z  ]  (2.10) 

— ;  n  n  n  n  n  n 

z  n~l 
n 

where  M  and  N  are  positive  definite  matrices  and  L  is  a  positive  semi  definite 

T 

matrix  such  that  for  C  C  =  L  the  pair  (F(p),0)  is  observable.  This  criterion 

gives  2  in  the  form  (2.8)  with 
n 

s  =  s  ( p )  =  [D(p)T.-(p)D(p)  +  N]_1  [D(PH(p)F(p)+M] 
where  :(p)  is  obtained  from  the  algebraic  Ricatii  equation 
S>(p)  =  K  (p)  +  L  -  ^2 (p)K^ (pjK^Cp) 

where 

K  ^ (p )  =  FT(p)$(p)F(p) 

K-2fp>  =  F(p)<t>(p)D(p)+M 
K3(p)  =  D  (p )  <£  (p  )D  (p  )+N 


/ 
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Example  1 


Consider  the  system  (2.7)  with 


F  (p)  = 


ri 


,  D(p)  = 


,  G(p)  =  I,  H(p)  *  [1  0] 


and  state  feedback. 


Z  =  -  Is.  sjx 
n  1  2  n 


so  that 


F(p,s)  = 


prsi  p2_s2 


It  can  be  seen  that  the  system  (F (p , s ) ,C(p) ,H(p) )  is  controllable  and  observable 
and  the  pair  (F(p),  D(p))  is  stabil izabley  independently  of  the  values  of  P1*P2*S1*S 
The  conditions  of  theorem  1  are  thus  satisfied  and  the  cost  function  J(p,s)  has  a 
finite  limit  value  for  all  p  and  s.  On  the  other  hand,  taking  F(p,s)  as  above  but 


/ 
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wi  th 


r 1 

1 

D(p)  =  ! 

,  G(p)  = 

1  0 

0 

H(p)  =  [0  1] 


it  can  he  seen  that  the  system  (F (p , s) ,G(p) ,H(p) )  is  controllable  and  observable 
everywhere  but  on  tiie  line  =  s^  and  the  pair  (F (p ,s) ,D (p) )  is  controllable 
(hence,  stabll  izable)  everywhere  but  at  p^  -  0.  Since  these  points  in  the  P  S 
hype rp lane  produce  indefinite  values  for  the  cost  function  J(p,s),  they  may  simply 

be  avoided  in  the  search  procedure. 


ft  should  be  emphasized  that  while  the  conditions  of  theorem  1  ensure  that  the 
fixed  gain  design  criterion  is  well  defined,  they  do  not  guarantee  that  the  resulting 
fixed  gain  controller  will  be  stabilizing  at  all  points  in  the  given  parameter  set  P. 
In  tact,  overall  stability  may  not  be  achievable  by  a  fixed  gain  controller  for  given 
parameter  sets.  To  ensure  stability  the  parameter  set  must  be  restricted  in  such  a 
way  that  the  eigenvalues  of  F(p,s)  are  mapped  into  the  unit  circle  for  all  values  of 
p  and  s.  This  problem  has  been  recently  considered  in  [6], 

C .  Filter  Design. 

Once  a  fixed  gain  controller  has  been  selected,  the  system  is  given  by  (2.1) 

o 

with  s=s  .  To  each  parameter  p  there  corresponds  a  linear  stochastic  model  and 
an  optimal  estimator  (Kalman  filter)  based  on  this  model.  It  is  now  desired  to 
find  an  estimator,  independent  of  p(i.e.,  non-adapt ive) ,  that  will  be  optimal  in 


T' 
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some  sense  for  the  entire  parameter  set  P.  Restricting  the  search  to  the 
set  of  Kalman  filters  defined  on  P,  the  minmax  information  criterion  defined  in 
Section  A  implies  that  the  desired  estimator  is  the  Kalman  filter  corresponding 
to  a  parameter  point  te.P,  which  is  found  from 

f  (Y) 

min  max  E  log  — ^ - 

teP  peP  P  f  (Y) 


where  f  (Y),  f  (Y)  and  E  are  the  probability  densities  and  the  expectation 

p  t  p 

associated  with  the  corresponding  parameters.  (note  that  if  t  is  identified 
as  the  filter  design  parameter,  then,  using  the  notation  of  section  A,  we  have 

ftrp)  =  fP  =  fP,  fC  =  fC  -  f  and  Et(p)  -  EP  =  E  ). 

P  P  Pit  p  pp 


For  an  observation  sequence  Y 


.N 


(yr 


.  .,y^)  let  us  denote 


N 

2  fn(Y  > 

V"-0  -s  Ep  loetiA 


then  the  filter  design  criterion  in  finite  time  may  be  written  as 


min  max  I  (p,t) 
teP  pcP  N 


The  calculation  of  IN(p,t)  is  similar  to  that  of  IN(p,s)  in  the  controller  design 
problem,  replacing  in  the  corresponding  equations  F(p,s(p)),  F(p,s),  K(p,s) ,B(p,s(p>) 
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and  B(p,s(p))  by  F(t,s°),  F(p,s°),  K(p,s°),  B(t,s°)  and  B(p,s°)  respectively, 
where  s°  is  the  selected  controller  gain  vector,  and  replacing  H(p)  by  [ H ( t )  -  l'(p)l. 

As  in  the  case  of  controller  design,  steady-state  estimator  design  offers  a 
computational  advantage  over  finite  time  design.  We  have  the  following  result. 

Theorem  2 

Suppose  that  for  each  peP  the  system  (F (p , s°) ,G(p)  ,  H(p))  is  stable,  controllable 

and  observable.  Then  the  cost  function  I  (p,t)  has  a  finite  limit  value 

N 

Proof :  Substituting  the  corresponding  matrices,  the  proof  is  identical  to  that  of 

Theorem  1 . 

It  should  be  noted  that  an  intermediate  result  of  the  proof  is  that  under  the 
condition  of  the  theorem  the  state  estimation  error 


l>  „(t,p)  =  E  {[x  -x  |  (p)][x  -x  I  (p)  ]  } 

n  t  n  n|n-l  n  n|n-l 


has  a  finite  limit  value  for  any  t  and  p.  This  means  that  although  the  non-adaptive 
filter  is  not  optimal  in  the  least-squares  sense  at  any  point  but  t,  its'  mean-square 
error  has  a  finite  limit  value.  The  fixed-gain  steady  state  estimator  is  obtained 


I 

f 
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f  rora 


min  max  I(p,t) 
1 1: P  pcP 


where  I(p,t)  is  the  limit  value  of  I  (p,t). 
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3 .  Fixed  Gain  Controller  and  Filter  Design  for  Aircraft 

The  proposed  design  procedure  has  been  applied  to  the  problem  of  finding 
fixed  gain  values  for  the  longitudinal  control  system  of  a  given  fighter  aircraft 
(the  YF-16  CCV  data  has  been  used  in  the  design).  The  linearized  short  period 
equations  are  given  by  (e.g.  17]). 


where  q  is  the  pitch  rate,  a  is  the  angle  of  attack,  w  is  a  wind  gust  disturbance, 

<5  is  the  elevator  angle  (no  actuator  dynamics  are  introduced  here,  for  simplicity) 
e 

and  C  is  a  zero  mean  unity  covariance  white  noise  component  of  the  wind  disturbance. 
The  coefficients  d  (p)  and  g^(p)  are  functions  of  the  flight  condition 

parameters,  p  «  where  h  is  the  altitude  and  M  is  the  Mach  number.  Pitch 

rate  and  normal  acceleration  measurements  are  used  for  feedback  purposes.  The 
measurement  equation  is 


where  a 


is  the  normal  acceleration,  U  is  the  forward  velocity  and  v  is  a  zero  mean 
white  noise  process.  Typical  sensor  accuracy  and  wind  disturbance  data  can  be  found 
in  the  literature  (e.g.  [8]).  The  open  loop  dynamical  system  was  found  to  be 
unstable  for  most  of  the  given  flight  conditions. 

For  digital  control  purposes  the  equations  have  been  disretized  at  0.02  spconds 
sampling  intervals.  The  resulting  system  has  the  form 


Xn+1  =  F(p)xn  +  D(P)<5e  +  G(p)!;n 

n 


v  =  H(p)x  +  v 
’  n  n  n 


As  a  first  step  in  the  design  procedure  an  optimal  controller  is  designed  for 
each  of  the  flight  conditions.  While  any  desired  criterion  could  be  used  at  this 
stage,  we  have  chosen  for  illustration  purposes  the  linear  quadratic  criterion  [81 


min 


z 


n 


n=0 


nz 

'Se 

n 

+  n 

i  n 

2 

+  2 

+  2 

Snax 

a“ 

nz 

6e 

max 

max 


It  has  been  found  that  for  all  given  p  values  the  pair  (F(p),D(p)  is  stabilizable 
and  the  pair  (F(p),C)  is  observable  where 
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K 1 ight 
Cond  . 

No  . 

alt . 

(ft) 

Mach 

No. 

Dvn . 

Press.  „ 

(lb/ft  ) 

s 

■-'l 

s 

q 

S 

1 

S.L. 

0.6 

363.35 

-1.377 

-2.048 

:  i 

-1  .  5242  ; 

’  | 

•> 

15,000 

0.6 

205.07 

-1 .121 

-1  .101 

1 

-0 . 04 13  j  !  I 

3 

30,000 

0.6 

107 .85 

-0.9603 

-0.608 

-0.0177 

it 

4 

S.L. 

0.8 

646.32 

-1.3916 

-3.395 

-0.0077 

; 

1 

5 

15,000 

0.8 

364.57 

-1.197 

-1.943 

-0.0072 

h 

30,000 

0.8 

191.73 

-1.079 

-1.115 

-0.0027 

7 

S.L. 

0.9 

818.00 

-1.3175 

-4. 307 

-0.0030 

1 

1 

8 

15,000 

0.9 

461.41 

-1.1903 

-2.575 

! 

0.0051  j  j 

9 

10,000 

0.9 

242.66 

-1 .1057 

-1 .525 

0.0028  !  j 

10 

S.L. 

0.95 

911.41 

-1.2 

-5.127 

0.0261  1 

1  1 

1  3,000 

0.95 

514.10 

-1 .1086 

-2.935 

i 

0.0456 

12 

30,000 

0.95 

270. 37 

-1 .0575 

-1 .687 

0.0225 

13 

S.L. 

1.2 

1  ,454.22 

-1 .6637 

-6.782 

0.1931  i 

1 

14 

15,000 

1.2 

820.28 

-1 .4532 

-4.044 

i 

0.2061  l 

| 

1  5 

30,000 

1.2 

431.39 

-1.3345 

-2.344 

! 

0 . 1092  ; 

1. 

16 

30,000 

1.7 

865.77 

-1  .7801 

-2.698 

i 

0.1281 

17 

50,000 

1.7 

335.37 

-1.68  37 

-1  .  383 

0.0554 

Table 

1.  Optimal  (Linear  Quadratic) 

Fecdii-!  k  <-a ins 

tor  Liven  Flight 

Com!  i  t  ions  . 

i 

i 

! 

! 

L 


-20- 


.so  that  the  resulting  controller  (regulator! 

z  =  -s(p)x 
n  n 

is  stabilizing.  In  addition,  the  system  (F(p,s),  C(p)  ,  H(p)l  is  control  i  abU-  and 
observable  for  any  s  and  p  so  that  the  conditions  of  Theorems  !  and  2  arc  sat  isfied 

and  the  fixed-gain  design  criteria  are  well  defined. 

In  searching  for  a  fixed-gain  controller  we  have  restricted  ourselves  ic  t he 
finite  set  of  optimal  gain  values  s(p)  for  the  17  given  flight  conditions.  These 
are  given  in  Table  1,  The  fixed  gain  controller  design  criterion  now  becomes 

min  max  ll(j,i)  ;  i.j  ■  1,...,17} 

i  .1 

The  resulting  fixed  gain  controller  gains  are 


s  =  -1.4532 

q 


s  =  -4.0437 
a 


s  = 

w 


0.2061 


-21- 


(these  are  the  optimal  controller  Rains  corresponding  to  flight  condition  14). 

It  remains  to  find  a  flight  condition-independent  state  estimator.  Substituting 
the  fixed  controller  gains  into  the  dynamical  models  for  each  of  the  flight 
conditions,  we  get  17  "closed-loop"  models.  Employing  the  filter  design 
criterion  presented  in  section  2.2,  and  restricting  the  search  to  the  set  of  17 
optimal  Kalman  filters,  the  filter  corresponding  to  flight  condition  3  is  found 
to  be  the  overall  (non-adapt ive)  optimal  filter.  The  resulting  control  svstem  is 
given  schematically  in  figure  3.1 
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Figure  3.1  Resulting  Fixed-Cain  Control  Svstem 
for  Aircraft 
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Although  we  did  not  specifically  introduce  performance  requirements  into 

the  design  procedure,  it  is  interesting  to  examine  the  performance  of  the 

resulting  control  system  at  different  flight  conditions.  Loosely  speaking,  it 

might  be  expected  that  with  a  fixed  gain  controller  the  aircraft  responses  would  be 

faster  at  high  dynamic  pressures  than  at  low  pressures.  (The  dynamic  pressure  is 
1  2 

p .  «*  -r  p  U  where  o  is  the  atmospheric  air  density  and  C  is  the  forward  velocitv.; 
d  l  o  o 

For  high  performance  military  aircraft  it  is  necessary  that  even  at  low  dynamic 
pressures  the  control  system  maintains  sufficient  maneouvring  power  (i.e.  fast 
response)  which  is  acceptable  to  the  pilot. 


/ 
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The  aircraft  response  in  the  pitch  rate  anti  in  the  normal  acceleration 
to  a  4  disturbance  ("initial  condition")  in  the  angle  of  attack  is  shown  in 
Figure  3.2.  A  linearized  aircraft  model  has  been  used  in  the  simulation. 

The  flight  conditions  considered,  ordered  by  increasing  dynamic  pressure  are 
3,  1,  5  and  11.  It  is  seen  that,  as  expectet i,  the  response  sneed  genera] lv 
increases  with  the  dynamic  pressure  (Note,  however  that  this  dependence  is  not 
simple.  For  instance,  the  response  is  significantly  faster  at  flight  condition  3 
than  at  flight  condition  1,  although  the  corresponding  dynamic  pressures  are 
nearly  the  same.)  Note  that  an  acceptable  settling  time  of  2-3  sec.  is  obtained 
for  the  lower  dynamic  pressure  range  (e.g.  flight  conditions  1  ant'  3)  while  for  the 
higher  pressure  range  it  is  much  shorter  (1  sec). 


■  -  ....  >  .  .. 
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’ure  3.2 


Pitch  rate  (left)  and  normal  acceleration  (right)  responses  of  resulting 
fixed-gain  control  system  at  flight  conditions  3  (top),  l,  “S ,  11  (bottom) 
The  dynamic  pressure  is  lowest  at  f.c.3  and  highest  at  f.c.  11. 
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Next ,  the  fixed  gain  design  is  evaluated  using  a  pilots'  opinion  .hart 
proposed  in  [7],  The  damping  ratio  and  the  natural  frequency  of  the  closed- 
loop  system  at  each  of  the  flight  .  'nditions  have  been  calculated  and  placed  on 
the  chart  as  shown  in  Figure  3.  It  is  seen  that  the  system's  performance  at  most 
of  the  flight  conditions  is  rated  as  satisfactory  or  acceptable.  (The  tendency 
towards  the  fast  response  range  and  away  from  the  slow  response  range  is  also 
desirable  for  high  performance  fighter  aircraft).  This  is  particularly  remarkable 
in  view  of  the  large  variability  of  the  flight  conditions  and  the  corresponding 
aircraft  dynamics,  and  in  view  of  the  fact  that  performance  or  handling  qualities 
have  not  been  specifically  considered  in  the  design. 


h 


n.-mt  Ing  Ratio 


Figure  3.3  Performance  evaluation  of  proposed 
fixed-gain  design  at  given  flight 
conditions  bv  pilot  rating  contours. 


A  -  Satisfactory 

B  -  Acceptable 

C.  -  Poor 

p  -  Verv  slow  response,  large 

control  motion  to  manoeuvre, 
difficult  to  trim 

E  -  Initial  response  fast, 

oversensitive,  light  stick 
forces 

F  -  Sluggish,  large  stick  motion 
ant'  forces 

(I  -  Slow  initial  lv,  then 

oscillatory,  tendency  to 
overcontrol 
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C.oncluding  Remarks 

A  method  for  designing  fixed-gain  controllers  and  filters  for  stochastic 
systems  with  large  parameter  variation  has  been  presented  with  particular  attention 
to  the  control  problem  of  high  performance  aircraft.  The  resulting  design  has 
shown  good  performance  qualities,  although  performance  specifications  have  not 
been  introduced  directly  into  the  design.  This  may  be  done  bv  first  designing  each 
of  the  pointwise  optimal  controllers  (at  each  of  the  flight  conditions)  so  as  to 
meet  the  specifications,  and  then  imposing  the  specifications  as  constraints  on  the 
minmax  search  procedure.  In  this  paper  we  have  found  a  fixed-gain  control  svstem 
for  the  entire  flight  regime.  However,  since  the  fixed-gain  controller  would 
normally  he  used  as  a  back-up  system  in  emergency  situations,  a  small  subset  of  the 
flight  conditions  corresponding  to  those  situations  may  be  considered  in  the  design, 
so  as  to  meet  tighter  performance  requirements. 
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